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VIBRATIONS OF BEAMS ON MANY SUPPORTS 
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(Proc. Paper 863) 


SUMMARY 


The natural frequencies of a continuous beam resting on an arbitrary num- 
ber of uniformly spaced supports are determined from a difference equation 
formulation. These frequencies fall in periodically spaced groups that are 
separated by spectral gaps of widths equal to approximately half the interval 
between the natural frequencies of a single beam on a square root frequency 
scale. These groups tend to uniform spectra as the number of supports tends 
to infinity, but the gaps remain, giving a band-pass character to the entire 
spectrum. Wave propagation along an infinite, periodically supported beam 
is discussed and the phase and group velocities evaluated as functions of fre- 
quency. 


INTRODUCTION 


The natural frequencies of continuous beams on several, uniformly spaced 
supports have been determined by Ayre and Jacobsen,(1) their results being 
presented in the form on a nomograph. The purpose of the present paper is 
to present an analytical solution, which offers certain advantages, especially 
when the number of supports is large. The limiting case of an infinite beam 
on periodically spaced supports is of particular interest, both because of the 
band-pass character of its spectrum and because it has been suggested as a 
model for the study of flutter of periodically supported aircraft structures, (2) 
albeit without recognition of the nature of the resulting spectrum. 

We remark that both the problem and the method of solution considered 
here** find their natural antecedants in Brooks Taylor’s classical study (1713) 
of a massless, elastic string carrying equally spaced, point masses.(3) The 
Note: Discussion open until May 1, 1956. Paper 863 is part of the copyrighted Journal 

of the Engineering Mechanics Division of the American Society of Civil Engineers, 

Vol. 82 No. EM 1, January, 1956. 

* Prof. of Eng., Univ. of California, Los Angeles, Calif. 
** The application of finite difference methods to beam loading problems is, 
of course, well known; see, e.g., Th. v. Karman and M. Biot, Mathematical 

Methods in Engineering (McGraw-Hill, New York, 1940), p. 442. 
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analysis of such systems was further extended by Rayleigh(4,5) in order to 
generalize the wave propagation phenomena suggested by Reynolds’ discon- 
nected pendulums and the mechanical models of atomic radiation suggested 
by Fitzgerald (although Rayleigh himself pointed out that such models gener- 
ally were not consistent with the occurence of the first power of the fre- 
quence in spectral formulae such as those of Balmer). 


Formulation of Problem 


The equation of motion of a uniform beam is 


4 
(1) 


where 


E = modulus of elasticity 

I= moment of inertia of beam cross section with respect to neutral axis 
™ = mass per unit length of beam 

y = transverse displacement 
x =z 
t = 


coordinate measured along beam 
time 


Asmnming & harmonic time dependence of angular frequency W, we may re- 
place (0 y/d t2) by -W 2y. In the formulation of the boundary value problem 
for a beam supported at x = 0, (, 2 {, - - N ¢ (Fig. 1) we also find it con- 
venient to introduce a dimensionless, local coordinate such that 


y= yn(€), nel, (2a) 


€= ~(n-1), (n-l)2 x=nk (2b) 


Introducing, in addition, the dimensionless, eigenvalue parameter 


f/EI (3) 


(1) goes over to 


x 


x= O 20 (N-2)¢ (N-I)f 
2 N 


Fig. 1. A beam resting on uniformly spaced supports. 
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n=1,2,°°°N (4) 


We shall assume that the beam is simply supported at the end points as 
well as the intermediate points, although the following method of solution is 
equally amenable to other end conditions (the modification of the intermediate 
support conditions,* on the other hand, would be rather more involved). The 
boundary conditions then read 


Yn(O) = (5a) 
yn-(1)= y\(O) , (5b) 
, one (5c) 
(O)= (5d) 


Solution 


The most general solution to (4) that satisfies (5a) is given by 


= A, [sinh sin - sind sinh 
+B, [sinh d sin A(I-E)-sin sinh A(l-€)] (6a) 


In order to avoid indeterminate forms, we find it convenient to treat separate- 
ly those modes that are identical with those of a single span of length { - viz., 


y(€)=sinhk€, sind=O (6b) 


which clearly satisfy all of (5a,b,c,d). 
We first impose the moment continuity condition (5c) on (6a) to obtain 


B= AL, (7) 


Substituting (7) in (6a) and imposing the slope continuity condition (5b) on the 
result then yields the difference equation 


+A,_, =O (8) 


* E.g., to elastic supports. We remark that clamping the beam at the inter- 
mediate supports would render the motion of the successive bays independ- 
ent. 
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where we have introduced the parameter @ according to 


cos = (sinh X-sind) (9) 
The most general solution to (8) is given by 


A, = A sin(n@) + B cos(n@) (10) 


where A and B are arbitrary constants. The boundary conditions to be satis- 
fied by this solution follow from the imposition of (5d) on (6a) and (7) and read 


A,=A,20 (11) 


Setting A, = 0 in (10) yields B = 0, whence Ay = 0 requires sin(N @ ) to vanish. 
We are free to choose any value of @ that satisfies the latter condition, but in 
what follows we find it expedient to choose 


= (12a) 
(12b) 


where r is any positive integer, and (12b) implies @ = @ , in (9). 
The solution to the boundary value problem presented by (4) and (5) now 
may be placed in the form* 


)= sin(n§,) £(é )+ sin(n-1)4 (1-€) (13) 
£.(€) = sinh A, sind,€—-sin A, sinh (14) 


This form of the solution is indeterminate if A is an integral multiple of 77, 
but then we may resort directly to (6b). 


The Frequency Spectra 
We now proceed to examine those values of A determined by (9) and (12), 
which, in turn, determine the permissible frequencies of free vibration via 
(3). We first remark that 7 is the smallest value of A for which @, as given 


by (9), is real -- viz., cos@ = -1, corresponding to which we choose O=7 
and r= 1. Setting r = 1,2,3,—N then yields the successive values 


Ay <1.50567 (15a) 


* Illustrative plots of these mode shapes are given in footnote ref. 1. 
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Setting r = N + 1 gives a value of A between AN and 1.505677, but the corres- 
ponding solution given by (13) vanishes identically and is discarded as trivial. 
There are, however, no values of A in the interval (1.505677, 27) that yield 
real 6, and we find as the next sequence 


2m = < Anse (I5b) 
etc., the general sequence being given approximately by* 


ST = Ans + 2(- 


These conclusions are illustrated graphically in Fig. 2, the numerical values 
shown there corresponding to N = 6. 

The large values of the hyperbolic functions for A >7 naturally suggest the 
use of the asymptotic approximation (differing from the exact result (9) by at 
most 0.02) 


(I6a) 


or, equivalently, 


sin 2\ = sin’ 8, (I6b) 


In this approximation we have only to determine the first sequence Ay, ho, 
--Ay from (16b), after which the remaining A ,. satisfy the (approximate) 
periodicity condition 


Semi-infinite and Infinite Arrays 


The results of the preceding section remain valid as N tends to infinity 
(semi-infinite array), but then @ may assume any value, a and A may assume 
all values for which the magnitude of the right hand side of (9) is smaller than 
unity. The sequences of (15) then go over to blocks of continuous spectra - 
viz., 


* The right hand, upper bounds in (15c) follow from setting @ = (s + 1)7 in 
(9) and correspond to the natural frequencies of a single span beam with 
clamped ends. 
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We remark that the character of the spectra represented by (18) might have 
been anticipated by analogy with the properties of band pass filters. In con- 
trast to the spectrum of a filter made up of discrete elements, however, the 
number of pass bands for the continuous beam is infinite, although there can 
. no free oscillation below the minimum cutoff frequency corresponding to 
= e 
If we assume the array to be fully infinite (periodically supported, with 
N = -co, ---, -1,0,1, ---oc) the spectra are still given by (18), but we no longer 
are required to reject the solution cos(n@ ) to the difference equation (8) and 
must add to (13) a second, independent solution in which cos (n 6) replaces 
sin(n@). Alternatively, we may write 


ym = + (19) 


where f )(§) is given by (14) with A depending on @ through (9). 

Restoring the time factor ei“t, the solution (19) represents a disturbance 
travelling in the negative or positive x direction as @ is positive or negative, 
respectively. The wavelength of this disturbance is given by (6) 27/6 and is 
integrally related to { only when @ is an integral multipie of 7. It follows 
that the phase velocity is given by(7 


V= wh/8 (20) 


where, we recall, 8 depends on w through (3) and (9). The group velocity then 
is given by(8) 


= 


The derivative (dw/d@) may be evaluated from (3) and (9); it is, of course, 
real only when @ is real but then has the same sign as @ — i.e., the phase 
and group velocities are of the same sign (that U and V could be of opposite 
sign is not inconceivable). The corresponding values of U and V for an infi- 
nite, unsupported beam are given by [let y = ei (t-x/v) in (1)]_ 


Veo = = w? = (ELV (22) 


»/8 (23) 
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U/U,,= dd/d8 (24) 


On the basis of the approximation (16) these ratios become 
W/V, = d/sin" (sin (25) 


U/U,, = (sin (cos + sin A) (26) 


The approximations (25) and (26) are plotted in Fig. 3. We observe that the 
phase velocity ratio decreases monotonically from 1 to 1 - (2s + 2)-1 asA 
increases from s7 to (s + 1/2)7, whereas the group velocity ratio vanishes at 
both ends of these intervals and achieves the maximum value 2 -1/2 at the 
midpoints of the intervals. These properties are again generalizations of 
their counterparts for a discrete element, band pass filter; in particular, the 
vanishing of the group velocity is associated with the existence of standing 
waves corresponding to the independent vibrations of the individual spans at 
frequencies corresponding to pinned or clamped ends as A = sm or (s + 1/2)z, 
respectively. 

We conclude by noting that properties similar to those discussed in the 
present paper are characteristic of any continuous system subjected to inter- 
mediate constraints, a fact first noted by Kelvin (1881): ) Numerous examples 
are given by Brillouin(10) but the problems associated with beams are unique 
insofar as they are governed by a differential equation of fourth order. 


February 19, 1955 
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ON THE DEFLECTIONS OF BOW GIRDERS OF 
NON-CIRCULAR SHAPES 


Enrico Volterra M. ASCE 
(Proc. Paper 870) 


SYNOPSIS 


Deflections of bow girders of non-circular shapes built-in at both ends 
are studied. Solutions are given in explicit form in the cases in which the 
central line of the curved beam is a cycloid, a catenary or a parabola and the 
external forces are uniformly distributed along the whole length of the beam. 
Curves of the deflections, angles of twist, bending and twisting moments are 
presented for certain cases and are compared with the corresponding values 
for the circular bow-girder. 


INTRODUCTION 


Notation—The letter symbols in this paper are defined where they first 
appear, in the text or by diagram, and are assembled alphabetically in 
Appendix I. 

In Chapter Nine of his recent book entitled “Studies in Elastic Struc- 
tures”2 (1), Professor A. J. S. Pippard discusses many interesting cases of 
circular arc bow girders under distributed or concentrated external loads. 
His method of approach to the problem is as follows: By expressing the total 
strain energy due to bending and twisting only, since the contribution due to 
shear can be neglected, and applying the principle of minimum strain energy, 
the three redundant unknowns, i.e., the bending and twisting moments and 
shearing force in the middle section of the arc, are calculated. 

If the external load is divided into a symmetrical and a skew-symmetrical 
system of loadings, the application of the method of superposition simplifies 
the numerical computations involved in solving the three simultaneous 
equations. 

In fact, only one equation has to be solved for the bending moment, and 
two simultaneous equations for the twisting moment and shearing force. 


Note: Discussion open until May 1, 1956. Paper 870 is part of the copyrighted Journal 
of the Engineering Mechanics Division of the American Society of Civil Engineers, 
Vol. 82 No. EM 1, January, 1956. 

1. Professor of Mechanics, Rensselaer Polytechnic Institute, Troy, New York. 


2. Numbers in parentheses refer to the Bibliography at the end of the paper. 
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Summing up algebraically the results of the two cases dealt with separately, 
the final required solutions are obtained. 

But in the more general cases of non-circular arc bow girders fixed at 
both ends Professor Pippard finds that the resulting expressions “cannot 
usually be integrated by exact mathematical process but must be calculated 
approximately by graphical or arithmetical methods.” 

However, it can be shown that also in the case of a non-circular bow 
girder built-in at both ends and subjected to a uniformly distributed load the 
results can be expressed in closed form if the central line of the curved 
beam is a cycloid, a catenary or a parabola. 

In order to obtain these results, a different and more general method of 
approach than the one described by Professor Pippard in his book, has to be 
applied. In fact, four differential equations connecting the bending moment 
Mx, the twisting moment Mz, the deflection v, the twisting angle f, the initial 
curvature R of the central line of the beam, and its inertia characteristics 
El, and torsional characteristics K have to be used. 

Two of these equations are the well-known De-Saint Venant’s equations 
for a curved beam bent out of its plane of initial curvature(2,3), while the 
two subsidiary equations of equilibrium are obtained by expressing the con- 
ditions of equilibrium of the beam. 

The method just described is the same as the one already used in previous 
papers in the study of problems of deflections of circular rings on rigid sup- 
ports and circular arc bow girders under distributed or concentrated 
forces(4), and of circular rings supported by elastic foundations and subjected 
to concentrated symmetric or antisymmetric forces(5,6,7), 

But, due to the fact that in the present case the curvature R of the central 
line of the bar is not constant, the solution of the system of differential equa- 
tions is more laborious than in problems of circular beams. However, since 
in the particular cases in which the central line of the beam is a cycloid, a 
catenary or a parabola the relationships between radius of curvature R ina 
generic point of the central line whose normal makes the angle w with the 
axis of symmetry of the beam (see Fig. 1) and the curvature Ro at the central 
point P, are expressed by the following simple formulas: 


for the cycloid: R= cos 


R 
for the catenary: R = —2 
cos 
for the parabola: R= 
cos* 


and if the external forces are uniformly distributed along the whole length of 
the beam, then the results for the bending and twisting moments, deflections 
and angles of twist can be expressed in closed form. 

Graphs of the deflections, angles of twist, bending and twisting moments 
are presented for certain cases and are compared with the corresponding 
values for the bow girder of circular shape. 

In the preparation of these graphs the author acknowledges with thanks 
the help given by Messers T. E. Edelbaum and E. C. Zachmanoglou, both 
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students in the Department of Aeronautical Engineering at Rensselaer 
Polytechnic Institute, and the help by the “Istituto per le Applicazioni del 
Calcolo” of the Italian “Consiglio Nazionale delle Ricerche,” Rome, Italy, in 
connection with some of the numerical computation involved. 


The Equations of the Problem 


The beam will be referred to a system of rectangular coordinates OXYZ 
with the origin O at the centroid of the cross-section, the axes X and Y 
coinciding with the principal axes of inertia of the cross-section, and the 
axis Z with the tangent to the center line at O, as shown in Fig. 2a. If My 
and Mz are the moments acting on the cross-section at O about the X and Z 
axes respectively (bending and twisting moments on the generic section of 
the beam), v the displacement of the centroid O in the direction of the Y axis, 
# the angle of twist of the cross-section about the Z axis, El the flexural 
rigidity, K the torsional rigidity, R the initial radius of curvature of the cen- 
ter line of the bar, the following two equations will be verified(2,3): 


a 1 dv 
wo = k| 
as RdS 


The equations of equilibrium are: 


as *dS R 


where p represents the external load, which is supposed uniformly distributed 
along the whole length of the bar.* 


*Equations (2) can be obtained by expressing the conditions of equilibrium of 
an element of the bar of length dS. In fact the following equations can easily 
be derived from Figure 2b: 


da 
2 sin —+ - +d! 0 
but, since: 


2 sin 2% = aa = 9S 
2 R 


M,.da - dM, = 0 


3 
7 
= 
as* 
(1) 
= 
v4 
Aig (2) 
ds R 
_ 
|: 
16 
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In the case in which the arc of length 2 § is symmetrical with respect to 
point P, (see Figure 1) the conditions of fixity of both extremities of the 
beam imply the following boundary conditions: 


ds 


= viz2) = (3) 


Introducing as new variables the adimensional quantities: 


PRo 


*Continued from page 3. 


from which the equation Mx - uz = 0 follows. 


V (V+4dV) + pdS = 0 


+ VdS - M ein da = 0 


- aM, + VdS - M, da = 0 
aM, 
R 


as Zas as 


M 
r 
2 
x=8 (4) 
y= i dv 
ds 
where: 
A= = (= +—) 
2 k EI, 
av 
ds 
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and assuming as independent variable the angle yw that the normal to the 
central line of the beam in a generic point P makes with the axis of sym- 
metry at P, (see Figure 1), equations (1) and (2) become: 


+ - = 0 


ay 


*Continued from page 4. 
Equations (2) could also be obtained by minimizing the potential energy of 


the bar, which, considering only the terms due to bending and twisting mo- 
ments, since, as pointed out by Professor Pippard the term due to shear 


can be neglected, will be expressed by: 


5 
“se 
(8) 
‘ 
de £ 
where: 
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and C is a constant of integration to be determined from the boundary 
equations (3) or: 


B(+a) = v(t+a) = y(+a) = 0 (6) 


Case of the Circle 


In this particular case: 
p=1 


and the solution of equations (5) and (6) gives: 
m= 1+ A cos 


+A siny? 


2 
X= [2 + Aypsiny- cos y ( + Aatana 


A { - (u'+1) p [ | cos +2c0sa } 


p! = 


A = sine 


in accordance with the results already obtained in a previous paper (4) but 
using a different notation. 


Case of the Cycloid 


In this particular case: 


cos 


o = siny 


and the solution of equations (5) and (6) gives: 


1 1 
= cos A cos 
™ 3 4 


| 
et 
1 
(14u") 
. 
where: 
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m, =~ doin 2y+ easing 


x= 2. { cos sy - (34p) A cos 29+B cosy = (3-p) Ysiny + 
2 


+ 3(1-p) a} 


= 49 + 8Y+ 8 cosy) - 


-[s (1-8n)] cos 2+ sin 29) + 


B = pA(cos3a+3cosa) (38-cos4a) (1+p) cos2a 
cs re cos A(3=p)(cos8a+8cosa) + 
5 
+f (1-94) | cos2a (3-,)a(a+sin2a) 


Case of the Catenary 


In this particular case: 


and the solution of equations (5) and (6) gives: 


1 1 1 l+einy 
= Aco - +=- = siny log 
™ 2008 f 1-sing 


m, = A sing+ 2 tan y+ 2 cos log 
2 1-siny 


Ese 
7 
f 
: 
A = a 
73, 
where: 
1 
2 
cos*y 
o = tan 
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x= sing+ cos log|cos g| + B cos P+ 


2-3 2-3 i 
+ siny log cos? + ) cing ] 


u =~ loglcosy! _ yiog g + 
cosy cosy 2 l-siny cosy 


log? itsing 1-p sunt log i+siny 
240084 32 l-siny 8 cos 1-siny 


1 1 i+sin ] 
where: 


1 sina + 42 tog 1+s8in 2+» F,(a) 


A = 4 1-sina 
2ua- (1+p)tana 


B=-afi+p-2 106 | | « 


tandlog Atsind _ F, 
1-sind 


+ 


, logicosAl _ 1 Atsind , F,@) | 


cosa cosa 2 l-sind 


+ 5-7 , log? sind log _ 
a, 8cos%y 32 l-sina 8 cos% l-sind 


‘ - 4 F,(%) 10g + 3(4) | 
cosa 2 1-sind 


1-sinp 


as 
i 
: 
ab 
Rigs 
B 
C= 
1 l+sin 
FL“) = 4 log dp 
2 
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= 2 | tan p log A+sinP 4p (7) 
2 1-sinp 
=i tan p log* i+sinp dp 
4 1-sinp 


Case of the Parabola 


In this particular case: 


P 


cos 2 1-siny 


and the solution of equations (5) and (6) gives: 


Scosy 24cos 


x = (1-3p) Ac siny + in. + B cosy o* cosy + 
cos 


24 cosy 120cos°y 45 cosy 


cos*¥ any]+ 2c08 *p 8 cos*y 


5 (3-19) sin ¥ 5(3-19p) o Sin 


+ 
128 48 cos® 192 cos*y 
, 208-859 , 2(1-3p) ¢ 
128c08 » 5760cos 45cos4y 


a 
j 
12 
q 
\ 
ate 
2 
> 
ip 
a 
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where: 
g(a), sina sina + tana[ 1 + 2] 
} Scos*a 4 cosa __60 coats cos 2a 


(1-3 )o(a)+2p 
cos 


a 


,, o(a) |? + ~Sing_ | ofa) 
cosa 4 8cos*a L cos*a 


24cos°a 12cos*a 45 


3 o(a) 
= (1-3p) af ova) tan a - cos *a 


+ + o(a) + 5 (3-19 o(a) sina 


cos~a 192 cos "a 


359-203 + 2(3p-1) 
128c08%a 5760cos°a 45cos*a 


In figures 3 to 10 the quantities myx, mz, x and u are given for the following 
values of the parameters = 30°, 60°; w= - 1/7, 4 =0, In figure 11 
the quantities mx, mz, x and u are given in the case of the cycloid and the 
circle for a = 90°; u=-1/7, H=0, w= 1/3. In figure 12 the functions 
are given for the parabola, catenary and cycloid, while in figure 13 the values 
of the functions F1(y), F2(y) and F3(y) are given. 


APPENDIX I. NOTATION 


The following symbols, adopted for use in the paper, and for the guidance 
of discussers, conform essentially with American Standard Letter Symbols 
for Structural Analysis (ASA-Z10. 8-1942), prepared by a Committee of the 
American Standards Association, with Society representation and approved 
by the Association in 1942. 


A = a constant 
B = a constant 
Cc = a constant 
E = the modulus of elasticity 


Fj, Fg, Fg = functions defined by equation (7) 


= 
: 
| 
te 
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1, = the rectangular moment of inertia 
= the torsional rigidity 
= bending moment 


= twisting moment 
_ Mz 
PRG 
= point of symmetry of the bar (defined in Fig. 1) 
= intensity of distributed load 
= the initial curvature of the center line of the beam 
= the initial curvature of the center line of the beam at point P, 
= arc length 


= displacement along y axis 
= rectangular coordinate axes 
= an angle defined by equation (6) and by Fig. 1 
= the angle of twist 
= the angle of twist per unit length 
1 
Lk * 


+ 
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Discussion of 
“A NEW APPROACH TO TURBULENT 
BOUNDARY LAYER PROBLEMS” 


by Donald Ross 
(Proc. Paper 604) 


DONALD ROss.!—In attacking the problem of turbulent boundary layers, 
and in writing the subject paper, the writer has been primarily concerned 
with the case of boundary layers acted upon by adverse pressure gradients. 
For this reason, and also because several papers in the literature handle 
the subject successfully, that part of the original thesis dealing with equi- 
librium flows was omitted from the paper. However, in three discussions 
(see Proceedings Separate 776), Messrs. Baines, Robertson and Bauer have 
dwelled on this case and other readers may feel that it should receive 
specific attention. 

For equilibrium boundary layers, there seems to be little question that 
equation 5 of the paper correctly represents the velocity distribution of the 
inner turbulent flow, although there may be some variation of the coefficients 
in equation 6. However, the velocity deficiency in the outer region may more 
properly be represented by the more general equation B proposed by Baines. 
The writer agrees with Baines and Bauer that a higher exponent is in better 
agreement with flat-plate data than the 3/2-power used by the writer. The 
writer chose the three-halves power because for flows having strong adverse 
pressure gradients the three-halves power law is a better fit than any other. 
Also, because the velocity deficiencies involved are so much greater near 
separation, the writer settled on the single formula giving the best fit for 
that type of flow and accepted the error typified by Bauer’s Figure 14. It is 
not a serious discrepancy. In Figure 15 we may compare the writer’s 
formula: 


1-- =84Y— 
a 8.4 5 3 
with formulas of Baines and of Hama” and with the expression derived by 
von Karman that the writer finds gives the best experimental fit. This 


1. The Bell Telephone Laboratories, Inc., Whippany, N. J. 

2. Hama, Francis R., “Boundary Layer Characteristics for Smooth and 
Rough Surfaces,” Transactions, Society of Naval Architects and Marine 
Engineers, Vol. 62, 1954. 
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figure clearly shows that the differences between the writer’s formula and 
Baines’ are not great and that both give good average fits to experimental 
data. However, the true curve is not a power law, hence Bauer’s result 
(Figure 14) that the value of D is not constant across the profile. Baines 

has incorporated the coefficient cg in his formula for the outer velocity 
deficiency. This is valid for the constant pressure case, and the author 
confirms that D is proportional to + for flat plates. However, in 
general, for non-equilibrium boundary layers, i.e., boundary layers having 
profiles that change with distance, the local value of D is in no way dependent 
on the local value of Cr. 

In reply to Baines, the writer believes that the methods developed in the 
paper can be extended to rough surfaces to the extent that these are under- 
stood. Thus, following Hama“ and Clauser3, we may write for the 
logarithmic inner-turbulent velocity profile: 
= Blog 


u, 


AU 
where: 

AU 


u, 


= Blog za, + const. 
v 


in the fully developed roughness region, and is somewhat different in the 
transitional region. This factor is directly related to the increase of the 
surface friction that the roughness produces and its determination is dis- 
cussed by Hama.” Fora rough surface, equations 18 through 20 are ap- 
plicable if u Oe is repiaced by (ug, + AU). With this extension, most of 


the rest of the paper should be directly applicable to rough as well as smooth 
surfaces. As noted by Bauer, D will generally be larger for flows having 
rough walls. 

Both Messrs. Baines and Robertson have pointed out that 6/@Disa 
function of c¢ and/or Rg for flat-plate conditions, In equation 15 of the 
paper, the relationship is only intended to be approximate and will be so 
corrected in the Transactions. 

Another misprint has led Robertson to an erroneous analysis. Equation 20 


should read: 
uw u 
i: = 0.7 + 5.0 log Rg be 
Cf ug uy 


without the plus sign that appears in the original printed version. Robertson’s 
friction analysis based on that version must be in error and his simple con- 
clusion does not hold. 

Robertson’s discussion of separation, on the other hand, should be valid, 
at least to the approximations made for his final equation. As Ax can be as 
great as 500 times §j;, and as A can actually be much greater than he as- 
sumes, the term involving these two factors should not be dropped relative to 


{;. In fact, the Sidney University experiments of Newman mentioned in his 


3. Clauser, Francis H., “Turbulent Boundary Layers in Adverse Pressure 
Gradients,” Journal of the Aeronautical Sciences, Vol. 21, 1954, p. 91. 


4 
we 
| 
é 
2 
3 
ok 
: 
x 
q 4 


Paper 876 DISCUSSION 3 


discussion are an example where this term must be included. The curves 
shown in Robertson’s Figure D are the limiting values for rapid decelera- 
tions. 

In response to a request by Robertson, the accompanying table lists the 
sources of the experimental data and identifies the symbols used in Figures 9 
and 11, 

The writer wishes to thank the three discussers for their constructive 
comments, and he hopes that this closure sheds some light on the interesting 
questions that they. have raised. 


Aa 


7 8 9 
Figure 15 - Comparison of Outer Velocity Formlas for Flat Plates 
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Discussion of 
“FAILURE OF PLAIN CONCRETE UNDER COMBINED STRESSES” 


by Boris Bresler and Karl S. Pister 
(Proc. Paper 674) 


HENRY J. CowAN.!—The authors deserve to be congratulated on having 
reduced the problem of the failure of concrete under combined stresses to a 
criterion similar to that which has been found to govern the failure of ductile 
materials. This is a significant advance in the theory of the strength of 
materials. However, the authors dismiss too lightly some of the previous 
theories, notably Rankine’s maximum stress theory, and Mohr’s theory, of 
which Coulomb’s internal friction theory is a special case. 

While the mechanics of the failure of concrete are not fully known, it ap- 
pears to be always initiated by the formation of very fine tension cracks, 
which are invisible or barely visible, but can, by careful measurement, be 
detected.(18, 19, 20) The ultimate failure of the concrete may, however, 
take two forms: a cleavage failure, which results in a clean fracture and 
is associated with the existence of high tensile stresses, and a crushing 
failure, which is accompanied by an appreciable formation of debris and is 
associated with the existence of high compressive stresses. This ultimate 
failure is immediately obvious and leads to the collapse of a plain concrete 
specimen, either instantly or after a short period of time has been allowed 
for the completion of the process of failure. In a reinforced concrete 
specimen a cleavage failure can be prevented by inserting suitable tension 
or shear reinforcement; however, the strength of the specimen thereafter is 
dependent on the strength of the steel. 

It does not seem likely that a physical explanation of the process of 
failure in terms of atomic structure will be derived in the foreseeable 
future, since the structure of concrete, or even of cement, is vastly more 
complicated than that of a metal. A phenomenological explanation is neces- 
sarily only approximate, and it is desirable that it should be associated with 
clearly observed phenomena, Since in concrete there are two types of failure, 
it seems reasonable to use two criteria. The cleavage fracture is evidently 
dependent largely on the tensile strength of the cement, The aggregate plays, 
however, a great part in resisting compression, and the cement does not ap- 
parently alter its angle of internal friction appreciably at low triaxial pres- 
sures;(3,4)2 the failure of concrete in compression is, therefore, not unlike 
that of a highly cohesive granular soil. In putting forward this physical 
explanation, the writer is conscious that it is an over-simplification, which 
will give way to a more precise interpretation when the structure of concrete 
is better understood, 

The above considerations led the writer (21) to adopt the dual criterion of 
failure referred to by the authors. The maximum stress theory has been 
used on numerous occasions to account for the cleavage, or tension, fracture 


1. Prof. of Architectural Science, Univ. of Sydney, Australia, 
2. Authors’ reference, Separate No. 674. 
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of brittle materials, including concrete. Mohr’s theory, using a curved 
envelope, has similarly accounted for the crushing failure of concrete. The 
first attempt at a unified criterion of failure appears to have been made by 
Leon, 5) using Mohr’s theory with a continuous curved envelope, crossing 

the o - axis on the tension side. This is represented by the curved line ABA 
in fig. 8. Leon explained that a Mohr circle touching the envelope at B cor- 
responded to a cleavage failure, and Mohr circles touching the envelope else- 
where corresponded to compression failures. The curved envelope is, how- 
ever, an unsatisfactory criterion for the concrete designer, since it cannot be 
expressed in simple mathematical terms, and varies for different types of 
concrete. This difficulty is overcome by using a vertical limiting line DBD 
for the cleavage failure (Rankine’s maximum stress theory), and two limiting 
lines CC for the crushing failure. In the range of low triaxial stresses 
straight limiting lines may be used with sufficient accuracy (Coulomb’s in- 
ternal friction theory). Both the cleavage and the crushing failure can there- 
fore be expressed in terms of the principal stresses: a constant maximum 
tensile stress, and a maximum compressive stress which varies linearly 
with the minor (tensile of compressive) principal stress. The dual criterion 
of failure can be readily applied to reinforced and prestressed concrete de- 
sign. The linear relationship is sufficiently accurate for building construc- 
tion; it doee not apply to structures with high triaxial stresses, e.g., large 
concrete dams. 

Cowan and Armstrong (22)3 have checked this criterion in an investigation 
on reinforced and prestressed concrete beams under combined bending and 
torsion, and have suggested a procedure for structural design. The criterion 
has also been checked with the results of two earlier investigations on con- 
crete and reinforced concrete in combined bending and torsion, (23,24) In all 
cases, generally, good agreement with the experimental data was found. The 
only significant difference between theory and experiment occurred at the 
transition from primary torsion (cleavage) to primary bending (crushing) 
failure, since the dual criterion produces a sharp corner in the limiting line, 
whereas the experimental results follow a more gradual curve; however, 
this divergence is not of practical significance. 

A three dimensional representation in terms of the principal stresses of 
the criteria used by Bresler and Pister, Leon, and Cowan and Armstrong 
(Figs. 9 - 11) shows that the difference between them is due mainly to the 
effect of the intermediate principal stress: they are not as dissimilar as 
they appear to be when expressed in mathematical terms, and, indeed, this 
is borne out by a comparison with the experimental data. The intermediate 
principal stress does not play a significant part in a large range of combined 
stress problems, notably those met in the design of reinforced and pre- 
stressed concrete structures. Ease of application is, on the other hand, of 
prime importance in these problems; The writer therefore feels that the 
authors have been a little harsh in their discussion of some of the previous 
theories. This comment is not, however, intended to detract from the im- 
portant contribution which the authors have made to the theory of fracture of 
concrete. 


3. The full results of this investigation were published after the Authors 
wrote their paper. 
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Fig. 9—Bresler’s and Pister’s Failure Criterion. (0; and 
denote the principal stresses). (From Brown Boveri Review, 
August 1944) 


pein 


Fig. 10—Leon’s Failure Criterion. (o, 02, and 0. denote the 
principal stresses). (From Brown Boveri Review, 
August 1944) 
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BORIS BRESLER, ! A.M, ASCE, and KARL S&S. PISTER,?” J.M. ASCE.—The 
comments of the discussers fall into two main categories: (1) utility of pro- 
posed failure criteria for design, and (2) physical interpretation of failure 
criteria and mechanism of failure. Whether one form of failure criterion is 
superior to another is largely a matter of the use intended for it. The inter- 
action formula proposed by Professor Smith is useful in design but is en- 
tirely empirical, Professor Cowan proposes dual criteria which are some- 
what more rational but are more difficult of application, Neither criterion 
is applicable to triaxial states of stress. The writers’ criterion, although 
not directly applicable to design, may be useful for correlating biaxial and 
the triaxial tests and for furthering the understanding of the physical me- 
chanism of failure. It was the authors’ opinion that once a valid criterion 
is established, simplications suitable for design can be derived. For ex- 
ample, in the case of shearing stress T combined with compressive stress 
o equation (6) leads to the following definition of shearing strength of con- 
crete: 


The writers are grateful to Mr. Rice for pointing out the agreement be- 
tween the writers’ work and Professor Smith’s experiments. It should be 
noted, however, that the type of loading used in Professor Smith’s experi- 
ments produced a stress field which is difficult to analyze by elementary 
methods; therefore, the agreement obtained seems rather fortuitous. 

In formulating the failure law in terms of octahedral stress the writers 
were guided by several considerations: 


(1) The failure criterion in the form of Eq. (4), 7) = F(oo), accounts for 
the effect of the intermediate principal stress, G9, and thus is a natural 
generalization of the Mohr theory of failure which neglects the effect of ao. 
Furthermore, use of octahedral stresses allows plotting both biaxial and 
triaxial states of stress in one plane. A comparison of three possible forms 
of representation of a failure criterion is shown in Fig. 8. Figure 8a shows 
Mohr’s failure theory in terms of normal and shearing stresses; Figure 8b 
shows a failure criterion for brittle materials in terms of principal stresses; 
Fig. 8c shows a similar criterion in terms of octahedral stresses. Radial 
lines define some of the possible states of stress, and it is seen that only the 
octahedral form permits plotting triaxial states of stress in one plane. 


1 
= 0.10 + 8.9 (7) - 10.8 ey? 
c c 


(2) Octahedral stresses have physical interpretations which may lead to a 
better understanding of the mechanism of failure in concrete. The normal 
octahedral stress is the mean normal stress, and the shearing octahedral 
stress is proportional to the root mean shear stress at a point. Thus Eq. (4) 
suggests that the failure mechanism may be defined in terms of the mean 
shearing and normal stresses at a point. Another possible interpretation of 
octahedral stresses, as pointed out by Messrs. Blakey and Beresford, is re- 
lated to the elastic strain energy and the volumetric strain. It seems doubtful 
that elastic strain energy concepts based on homogeneous materials with rela- 
tively constant values of elastic modulus and Poisson’s ratio would be 


1. Associate Prof. of Civ. Eng., Univ. of California, Berkeley, Calif. 
2. Asst. Prof. oi Civ. Eng., Univ. of California, Berkeley, Calif. 
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applicable to concrete. This view is supported by the presence in Blakey and 
Beresford’s equations of odd terms, involving the second stress invariant Ip, 
which apparently have no simple physical interpretations. A valid adaption 
of strain energy criteria to concrete must be based on a physical concept of 
the mechanism of deformation and failure, 


(3) Within the range of states of stress considered by the writers, the 
octahedral stress criterion has a further advantage in its simplicity of form, 
resulting in a simple linear equation, the writers’ Eq. 6. 


Professor Cowan’s linear criterion of failure, differentiating between 
cleavage and crushing failures, is shown in Fig. 9 by dashed lines. Whether 
or not the transition from a pure cleavage to a shear fracture in concrete is 
gradual or abrupt seems to be a matter that is still unsettled. The parbolic 
envelope assumed by Leon for brittle materials provides a gradual transition 
from cleavage to shear fracture. For biaxial state of stress, 79 = 0, Equa- 
tion (6) may be rewritten in terms of principal stresses by substituting 
Eq. (5) into Eq. (6). This leads to a quadratic expression involving ¢o 

0 ¢, and the experimental constants. The equation may be cedaatntead bw ren rota- 
tion of the 04, 03 coordinate axes through a counter-clockwise angle of 45°, 
If the new coordinate axes are denoted by 0, 3 the following equation 
results: 
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This is the equation of an hyperbola in the pincipal stress plane, Fig. 9, and 
is similar to Fig. A by Messrs, Blakey and Beresford. It can be seen that 
Leon’s theory and the octahedral theory propose a single criterion applicable 
to both cleavage and shear failures, while Cowan proposes two distinct 
criteria with an abrupt transition between them. Aside from this aspect of 
their interpretation, the agreement between Cowan’s theory and that of the 
writers appears to be good. 

It must be understood that the proposed linear octahedral theory as yet is 
untested for states of stress in which the mean stress is tensile. (States of 
stress to the left of the 7, axis in Fig. 8c). It may well be that some modifi- 
cations in the assumed linear curve may be necessary. In order to test the 
validity of Eq. (6) beyond the range of the writers’ experimental data it may 
be of interest to extrapolate the straight line as shown in Fig. 10 and evaluate 
the strengths for various states of stress on this basis. The uniaxial tensile 
strength based on this assumption is 0.100, or 0.087 f'.. The torsional 
strength, or biaxial, equal tension-compression strength, is 0.100¢, i.e. very 
nearly equal to tension. The biaxial, equal tension-tension strength is 
0.07 o, or 0.061f'.., ie., seventy per cent of the uniaxial tensile strength. 
These values appear to be reasonable, 
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Messrs. Blakey and Beresford define failure as formation of microcracks 
at l. ads below the ultimate. The writers believe that the formation of micro- 
cracks is an important factor in the mechanism of failure, but that the extent 
to which these cracks jeopardize the load carrying capacity of concrete has 
not been definitely established. The work done by the discussers indicates 
ratios of cracking to an ultimate load varying from 0.43 to 0.80. Until the 
mechanism of failure in concrete is better understood the definition of 
failure shall remain arbitrary. 

The writers agree with Messrs. Blakey and Beresford that a complete 
study of strain distribution in the hollow cylinders would have been of 
interest. This study was not included by the writers primarily because 
interpretation of strain measurements in concrete is hazardous. Only sur- 
face strains can be measured without disturbing the internal structure of the 
material. Thus, the interpretation that deviation from linearity of load strain 
curve is evidence of microcracking is strictly applicable to the surface 
material only. The fact that surface cracking occurs at small fractions of 
the ultimate load seems to indicate that propagation of surface cracking with 
incrzase in load is slow, and that considerable load resisting capacity is 
developed after initial cracking. 

In conclusion the writers wish to convey their sincere appreciation for the 
valuable comments of the discussers, 
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Paper 876 


Corrections to 
“THICK RECTANGULAR PLATES ON AN ELASTIC FOUNDATION” 


by Daniel Frederick 
(Proc. Paper 818) 
Page 8, in Eq. 19, b, should be wp. 
Pp. 8 and 9, in Eqs. 19, 22, and 23, sin should be changed to cos. 
Page 9, in the first line beneath Eq. 23, 01,7 should be (2y. 
Page 10, at the end of the first sentence should be added the phrase: 


“and y and sin a» respectively.” 


Page 11, in the third line from the bottom, V,, (0) should be V, (0). 
Page 20, in Fig. 5d, C should be omitted. 
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